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Abstract Relying upon our previous treatment of the density matrices for nuclei (in general, nonrelativistic 
self-bound finite systems) we are studying a combined effect of center-of-mass motion and short-range 
nucleon-nucleon correlations on the nucleon density and momentum distributions in light nuclei ( 4 He 
and le O). Their intrinsic ground-state wave functions are constructed in the so-called fixed center-of- 
mass approximation, starting with mean- field Slater determinants modified by some correlator (e.g., after 
Jastrow or Villars). We develop the formalism based upon the Cartesian or boson representation, in which 
the coordinate and momentum operators are linear combinations of the creation and annihilation operators 
for oscillatory quanta in the three different space directions, and get the own "Tassie-Barker" factors for 
each distribution and point out other model-independent results. After this separation of the center-of-mass 
motion effects we propose additional analytic means in order to simplify the subsequent calculations (e.g., 
within the Jastrow approach or the unitary correlation operator method). The charge form factors, densities 
and momentum distributions of 4 He and ie O evaluated by using the well known cluster expansions are 
compared with data, our exact (numerical) results and microscopic calculations. 



1 Introduction 

Many efforts have been made to get a deeper understand- 
ing of the nuclear structure at small distances (less than 
the pion Compton wavelength) with realistic many-body 
calculations for the nuclear wave function (WF) whose 
short-range part strongly deviates from a mean-field de- 
scription. In this respect, as well known (see, e.g., survey 
[T] , ref . [S] and refs. therein) , the nucleon density matrices 
and their Fourier transforms are of great interest, being re- 
lated, on the one hand, to the nuclear ground-state (g.s.) 
properties and, on the other hand, to the cross sections 
of various medium- and high-energy scattering processes 
off nuclei. Regarding the second aspect, we mean firstly 
a comparatively simple relation in the Born approxima- 
tion to express the elastic electron scattering cross sec- 
tion through the charge form factor (FF) F c h(q) of the 
target-nucleus and its charge density p c h(f) being defined 
by the Fourier transform of F c h(q). In addition, in the 
so-called approximation of small interaction times (see [3J- 
|S]) the double differential (e, e') reaction cross section be- 
comes proportional to an integral of the momentum dis- 
tribution (MD) rj(p) over the momentum range that is 
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fixed with certain combination (the y— scalling variable) 
of the momentum transfer q and the energy transfer oj 
(cf. [6 ). Other links with rj(p) we find in approximate 
calculations of the spectral function that determines the 
exclusive A(e, e'N)X cross sections (see, e.g., review [7], 
ref. [S] and earlier papers [5], [TU]).Of course, two-body 
and more complicated reaction mechanisms, in particular, 
due to meson exchange currents (see, e.g., [IT] and |12| ) 
in electromagnetic interactions with nuclei, may obscure 
such links. 

Note also the distorted-wave-impulse-approximation 
calculations [T3] of proton MDs in 12 C and le O(e, e'p) 
reactions at Saclay kinematics, where the authors have 
shown a strong enhancement of the reaction cross sec- 
tions with account for the final-state interaction at re- 
coil momenta qn greater than 1.5 fm . In the range 
the corresponding distributions of outgoing protons, hav- 
ing a considerably slower fall-off with the g/{-increasing 
compared to the plane- wave-impulse-approximation ones, 
may imitate some SRC effect. Therefore, the correspond- 
ing theoretical approaches are needed in certain refine- 
ments to bringing a reliable information on the distribu- 
tions in question from experimental data. Neglecting these 
complexities one has to deal [3], [2], [IS] with the two 
structure quantities, viz., the intrinsic density distribu- 
tion (DD) or simply the intrinsic density pi n t(r) and the 
intrinsic MD rji nt (p). They are expectation values in the 
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translationally invariant (intrinsic) g.s. WF of appropri- 
ate many-body (multiplicative) operators which depend 
on the respective Jacobi variables. These definitions (see 
the next section) coincide with those by the Sapporo group 
|16| . [Tf] in studying the properties of few-body systems, 
but differ from the ones used by the authors of refs. [5], 
[TB] , [TS] in their calculations of the densities and momen- 
tum distributions in s — p and s — d shell nuclei. There 
we encounter the other (not intrinsic) quantities p(r) and 
n(k) introduced as in the case of infinite systems (e.g., 
the nuclear matter) by means of the expectation value of 
the one-body "density operator" with a trial Jastrow-type 
WF. The latter in its schematic form <P = F4> involves a 
correlation operator which incorporates correlations 
into the mean-field WF <P. It is required that F be trans- 
lationally invariant and symmetrical in particle permuta- 
tions. However, when starting with a Slater determinant 
(SD) <P, e.g., as in [2],[TB] the function & is translationally 
non-invariant ("bad"), that is, it contains spurious com- 
ponents which result from the CM motion (CMM) in a 
non-free state. In this connection, let us recall earlier and 
more recent attempts [20 -[27 to remedy such a deficiency 
of the nuclear WF, namely its lack of translational invari- 
ance (TI) wherever shell-model WFs (commonly built up 
from single-particle (s.p.) orbitals) are used. 

In most cases the CM correction has been made to 
calculate the FF F c h{q) and, respectively, the density 
Pchir) using, as a rule, the Tassie-Barker (TB) prescrip- 
tion (a comparison of the relevant effects can be found 
in ref. [15 ) while the not intrinsic DM n(k) has been cor- 
rected (without any good reasons) via the renormalization 

b — > ^j^b of the corresponding oscillator parameter b 

(see, e.g., [28]), i.e., as in the case of p c h{r). An alternative 
evaluation [3]-[3], [TS] of the intrinsic FF's, densities and 
momentum distributions, put forward in [3J to overcome 
some obstacles in describing the elastic and inclusive elec- 
tron scattering off the i He nucleus, has brought a fresh 
look at the CM correction of these quantities. In particu- 
lar, it turns out that Pi n t(r) and ri int (p) are shrunk (from 
the periphery of each of them to its central part) com- 
pared to p(r) and rj(p). To our knowledge, this significant 
consequence of the restoration of TI has been ignored in 
past and goes on to be missed in modern explorations [2], 

At this point, one should note that such a simulta- 
neous shrinking of the density and momentum distribu- 
tions has been found within the harmonic oscillator model 
(HOM) for the simple (Is) 4 configuration. Accordingly, 
the motivation of the present work is twofold. First, we 
will show our results obtained with WFs more realistic 
than the Is — shell SD composed of harmonic oscillator 
(HO) orbitals. Second, the approach of [5] ,[T5] is extended 
to heavier nuclei (cf. [14 ). The CM correction of their 
FFs and MDs is considered on an equal physical footing, 
viz., using one and the same translationally g.s. WF that 
incorporates the nucleon-nucleon short-range correlations 
(SRCs). 

1 Below, the notation F = C is employed as well 



We employ Jastrow WFs [29] and the unitary-model- 
operator (exp(iS) with = S) approach (UCOA) [3U] . 
[3T] and [32] to nuclear-structure physics and its develop- 
ment by the Darmstadt group [33J, [34 (cf. the diagram- 
free (coupled-channel) exp(S') - method with ^ —S 
in the many-fermion theory |35|V In the context, let us 
remind other methods of deriving the so-called cluster ex- 
pansions for the expectation values with respect to Jas- 
trow WFs (35]~ [13] • Among them we note a factor-cluster 
or Van Kampen-type expansion proposed in [41] to eval- 
uate the distributions of interest with special emphasis 
upon the correlated charge FF for elastic electron scat- 
tering off nuclei. It has turned out that the expansion 
is equivalent to an approximate version of the UCOA, 
described in |44] , and yields a factor-cluster analogue of 
the Iwamoto-Yamada expansion [36 . The former (called 
sometimes the FIY expansion) simplifies numerical calcu- 
lations compared to the latter. A careful comparison of the 
correlated one-body properties of s — p and s — d nuclei, 
evaluated within the Jastrow formalism by truncating the 
FIY, FAHT ( factor analogue of the expansion from |37| - 
[35]) and in the low-order approximation (LOA) from [32] 
for the one-body density matrix (1DM), has been carried 
out in |45| . In the three cases the CMM correction has 
been taken into account by the commonplace TB factor 
when extracting the model parameters (the HO parame- 
ter and correlation radius ) from the experimental charge 
FF (we will come back to the point later) . Of great inter- 
est are also the exact Jastrow calculations of the elastic 
FF, MD and two-body density of 4 He performed in [13] 
without any CMM correction (see our discussion below). 

The paper is organized as follows. The underlying for- 
malism with basic definitions is exposed in in the follow- 
ing section. Sect. 3 is devoted to constructing the trans- 
lationally invariant correlated WFs, while sect. 4 is con- 
tained the formulae obtained with the help of the UCOA 
decomposition of the similarity transformation C^O^'C 
truncated at the two-body terms. Here is a relevant 
one-body operator additive by nucleons. Explicit expres- 
sions for the DDs and MDs of nucleons in 4 He and le O 
are shown together with their FFs separately in subsect. 
4.1 and 4.2. Our results are discussed and compared with 
the data in sect. 5. Some intermediate derivations can be 
found in Appendices. 

2 The intrinsic form factor, density and 
momentum distributions and their 
evaluation in the Cartesian representation 

By definition, the intrinsic (elastic) FF of a nonrelativistic 
system with the mass number A and the total angular 
momentum equal to zero is 

1 A 

F(q) = F mt (q) = j^i&int | exp[zq • (r a - R)] | <F mt ) 

a=l 

(1) 

or 

F(q) = (& tn t | exp[iq • (n - R)] \ <F lnt ) = ... 



A. V. Shebeko et al.: Title Suppressed Due to Excessive Length 



3 



= (tfW | exp{iq ■ {r A - R)] | & int ), 

where Wi nt is the intrinsic WF of the system (nucleus), 
r a the coordinate operator for nucleon number a, and 
R = A- 1 J2t=i r a the CM operator. 

Recall that | Wi nt ) enters the eigenvector \\Pp) of 
the total Hamiltonian H of the system, which belongs 
to the eigenvalue P of the total momentum operator 
A 

P=EP a - 



one uses the terminology adopted in Chapter XI of the 
monograph [U]), which is related to the intrinsic MD |4J 



?7 int (p) EE A(\P int | fj int (p) | «P int ) 



(9) 



with 



\w P ) = \p) m, 



(2) 



Here p a is the momentum operator of the a-th particle. 
Henceforth the bracket | ) is used to represent a vector 
in the space of the center-of-mass coordinates, so that 
P\P) = P\P). A ket (bra) with an index | ••• ) a („(••• |) 
will refer to the state of the a— th particle. The intrinsic 
WF tf'int depends upon the A — 1 independent intrinsic 
variables. These may be expressed in terms of the Jacobi 
coordinates, e.g., 



%rt(p) = 5{P -Pa + P/ a ) = 5(p - flA-i) 

= \Va-i=p)(Va-i=p\- ( 10 ) 
The OBMD is the Fourier transform of the 1DM 

77 int (p) = (2n)~ 3 J d 3 rd 3 r'exp [ip • (r - r')]p£l(r, r'). 

(11) 



As in |15j we would like to point out that 
. 1 3 



Pint(r) 



A-l 



Pint \ 



"ii.tVA-l 7 *' A-l V ) 



(12) 



1 a 

Za=r a+ i--22 r P (a=l,2,...,A-l) (3) 

/9=1 

or the corresponding canonically conjugate momenta 

1 Q 
V a = — j(aP a +i - ^Pp) ( a = 1,2,..., A- 1). 

p=i 

(4) 

The WF Wpiri, r 2 , r A ) in the coordinate representa- 
tion satisfies the requirement of TI, 

V P (r 1 +a,r 2 +a, . . .,r A +a) = exp(iP-a)&p(n,r 2 , ...,r A ) 

for any arbitrary displacement a. 

The intrinsic density Pint(?~) is the Fourier transform 
of the elastic FF, or inversely, 

Fint(q) = jj e^V int (r)d 3 r. (6) 

From eq.Q it follows that pi nt (r) = A(&int\pint(r)\<Pint) , 
where 

A»t(r) = 6(r -r A + R) = 6(r - ^a-i)- (?) 
Further, the 1DM may be defined as 

p^ t (r,r')^A(^n t \p^(r,r'Mn t ) 

= A{* i n t \t A _ 1 = r}{t A _ 1 =r f \* int ) 

= Aj d 3 ^ . . . d 3 U-2#L (€i . • • • . r) 

X*int(tu---,tA-2,r'), (8) 

so that the normalization condition J d 3 rpH(7", r) = 
A is satisfied. We would like to emphasize that this is 
not an "imposed" definition. It appears naturally when 
evaluating the dynamical FF (or its diagonal part, if 



In other words, the intrinsic 1DM does not have the prop- 
erty pW(r) = pM(r,r) which can be justified for infinite 
systems, although it has often been exploited in approxi- 
mate treatments of finite systems (cf., however, ref. [48 , 
where an alternative definition of the 1DM for finite self- 
bound systems was proposed). 

Each of these intrinsic quantities can be written as the 
expectation value of a product of A operators acting on 
the subspaces of the separate A particles. For example, we 
have 

F mt (q) ee (& int | F mt (q)\V mt ) (13) 
with the multiplicative operator 



^int(<z) = exp[iqr • (f x - R)} = e 1 
whereas 



. . e 



Pint(r) = d(r 1 



R 



(27T) 



e-^ r F int (q)d 3 q. 



Now, we will use the Cartesian representation, in which 
the coordinate (momentum) operator r a (p a ) of the a-th 
particle is the linear combination of the Cartesian creation 
and annihilation operators a) and a . 



ro 



(fit 



P 



i^|(a f -a), r po = l, (14) 



with the Bose commutation rules 

[a], a}] = [a h a,j] = , 



[&i,a\] = 8, 



(15) 



The indices l,j — 1,2,3 label the three Cartesian axes 
x,y,z. 

As the "length parameter" ro one can choose the os- 
cillator parameter of a suitable HO basis in which the 
nuclear WF is expanded. Its basis vectors \n x n y n z )i ® 



\n x n y n z ) A , where the quantum numbers n x , n 



take on the values 0, 1 
\n x n y n z ) = [n x \ n y l n z l] 



are composed of the s.p. states 



|000), 
(16) 
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which are the eigenstates of the Hamiltonian H osc = 
cj(a + ■ d+ |), 

H osc \n x n y n z ) = (n x + n y + n z + §) u \n x n y n z ) , 

where uj is the oscillation frequency along the three axes 
x, y and z. We use the system of units with h = c = 1. 
The s.p. WF in coordinate representation is written 

(r | n x n y n z ) = Vv (x)ip ny [y)ipn s (z) , 
where (see, e.g., (49]) 

Ms) = [V^2"n!r ] _i H n {s/r ) exp(~s 2 /2r 2 ) 

and H n (x) is a Hermite polynomial. By definition, the 
oscillator parameter equals = [mu] ~~ 3 . 

Using eqs. (14- 15| , after some algebra one can get 



F lnt {q) = F TB (q) F H OM(q)x 

'A-r 



x exp 



x exp 



x exp 



x exp 



-n 



iq 



iq 

iq 

rp 

rp 



A 

A-l 
A 



rp -f 

7=2 a { 



V2A 



exp 



exp 



V2 



iq 



TP 

V2A 
ro 



a 2 



V2A 



(17) 



withF TB (g)= exp( 3 ^- g 2 rg), F HO M(q) = exp(-^q 2 r^). 

Thereat, the TB factor Ftb(q) appears automatically 
due to a specific structure of the operators involved. In 
other words, its appearance is independent of any nuclear 
properties (in general, properties of a finite system) . The 
only mathematical tool that has been used is the Baker- 
Hausdorff relation: 



,A+B 



e A e B e ' 



|[A,B] 



(18) 



that is valid with arbitrary operators A and B for which 



the commutator 



A,B 



commutes with each of them. 



2.1 Constructing intrinsic wave functions. Inclusion of 
nucleon-nucleon correlations 



A Slater determinant 
| Det) 



1 



e v V{\<j> Pl (l)) 

V€S A 



(19) 



as the total WF <P for an approximate and convenient 
description of the nuclear g.s., in the framework of the 
IPM or the Hartree-Fock(HF) approach exemplifies WF's 
which do not possess the property of TI, eq.([5]). Here e-p is 
the parity factor for the permutation V, 4> a the occupied 



orbital with the quantum numbers {a} and the summation 
runs over all permutations of the symmetric group Sa- 

There are different ways to restore TI if one starts with 
such a bad WF as | Det) ([50]-[52], (26|). 

According to Ernst, Shakin and Thaler (EST) pre- 
scription [5T] [^] in the fixed-CM approximation the nu- 
clear many-body WF with the total momentum P can be 
written in the form: 



EST\ 



| 9 P ) = \P) | * int 
The intrinsic WF after EST 



^est, _ (R = 0\$) 



[(<?|i?=0)(i?=0|#)] 1 /2 



(20) 



(21) 



is constructed from an arbitrary (in general, translation- 
ally non-invariant) WF <P, by requiring that the CM co- 
ordinate R be equal to zero. The corresponding FF is the 
ratio 

A(q) 



Fest(q) 



A(oy 



A(q) = (# | (27rr<y(Jl) exp[ i<7 • (fx - R)\ | (22) 
while the intrinsic MD 

, v (■P\(2nfS(R)6(p 1 -P/A-p)\-P) 

Vest(P) = 5 — ~ (2o) 

(<P | (2tt) 3 5(R) I 0) 

so that we have the Fourier transform 

Vest(p) = (2tt)- 3 J exp(-ipz)N(z)/N(0)dz (24) 
with 

N(z) = (<P | (2tt) 3 6(R) exp[z(pi - P/A)z] \ <P). (25) 

We see the certain resemblance between the structure 
functions N(z) and A(q), viz., both are determined by 
the expectation values of similar multiplicative operators 
with one and the same trial WF <P. Owing to this with the 
help of the same algebraic techniques (cf. eq. (17) ) we get 



A{q) = exp 



U(q) 



d\ exp 



with 



n 2 f 2 

q > o 



AA 
rp 



U(q), 
F(v,s) 



V2A 



(26) 
(27) 
(28) 



and the renormalized "length" parameter 



'•() 



A-l 



A 



-rp 



Other projection recipes can be applied without essential 
changes, see |15| 
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and. in parallel, Another popular option goes back to the lectures by 

Villars in jM] (see also [3T]) with a unitary operator 

^)=exp(-^K), (29) C = exp H G), (37) 

D(z) = [d\exp(- r ^r)F(v ',s% (30) £= (38) 

J \ 4A J a<0 

with where the Hermitian operator g(a,/3) acts onto the space 

p lr() of the pair (a,/3). In particular, we could follow the sim- 

z, v ' = (A - vp\z) (31) piest Darmstadt ansatz [33]: 



V2 ' V2A 



1 



and r- g(a,P) = (r a p)P a p+Pap* (*<*/»)}, (39) 

/A - 1 2 
Po = \/ — z — Po- 

V A where s is a function of the relative coordinate r Q( g = r a — 

When deriving these relations we have applied again rp. Its canonically conjugate momentum p aR = \(p a — 

eq.(18l in combination with the representation Pp)- 

Keeping in mind similar constructions we rewrite ex- 

(2tt) 3 s(r)= J exp (iXRj dX. (32) Pectation ^ as 

... ... . .. ... F(x,y) = ($(-x)\El(- V )E 1 (v)\$(x)), (40) 

Alter this we see that the expectations A(q) and N(z) are 

expressed through one and the same function F(x,y) where 

\$(x))=E 1 (x)...E A (x)\&), 
F(x,y) = ($\0 1 (x + y)0 2 (x)...0 A {x)\$), (33) „ . . 

since E x {x + y) = E x {x)E x (y) and \E a {x), Ep{y)\ =0 
where (a,j3 — 1, . . . , A) for any vectors x and y. 



1 {x) = exp(— £c*djj.) exp(a;a 7 ) = El(— x)E^(x) (34) 



Moreover, we find that 



E(x) r E- 1 (x)=r+ T ^ x (41) 
(7 = 1,..., A). V2 

In other words, we have constructed the generating func- an H 

tion for both. One should stress that this result has been E(x)p E^ 1 (x) = p — i^— x. (42) 

obtained independently of the model WF <P. 

1 dTbo? LTwF 11 PraCtiCe US C ° nSider a COTre " Remind tlmt Ef * In ° ther W ° rds ' Ea{x) iS the 

' displacement operator in the space of nucleon states with 

the label a. 



<£) =\$corr) = C '(1,2,- ••, A) \Det). (35) 



Due to this property when handling the similarity 
transformation 



x 



The ^4-particle operator C = C(r a — rp, p a — pp) ^intro- 
duces the SRCs and meets all necessary requirements of qi _ E x (x) E a{x)C{t — rp P ~ P 
the translational and Galileo invariance, the permutable 
and rotational symmetry, etc. However, being translation- xE~ 1 (x) E~ 1 (x) 
ally invariant itself such a model introduction of correla- 1 A 
tions does not enable to restore the TI violated with such we get 
a shell-model WF as the Slater determinant. 

What follows can be used with the Jastrow correlator O = C(E a (x)r a E~ 1 (x) — Ep(x)rpEp (x), 

m 

A J r A f/ - s , ols , E a {x)p a E- 1 {x)-E p {x)p p Ep 1 {x)) = 

C = -7g=> J=[[fir a p) (36) 

J a <P = C(r Q -r p , p a -pp) = C 

The normalization constant Cj = (Det | J* J | Det) (in i.e., 

general, a constant (Det \ C'G | Det), if any ) may C' = C. (43) 

be omitted keeping in mind the ratios A(q)/A(0) and 

N(z)/N(0). The function f(r a p) of the distance r a p = Reca11 that C is a Unction of all the relative coordinates 

\r a - rp\ is required to come to zero when particles a and and their canomcaUyconjugate momenta. 

P are inside a correlation volume of a radius r c . From e( l s - fl§5| and {43} it follows that 



3 



Of course, the operator may be spin and isospin dependent | <P corr (x)) = Ei(x) . . . Ea(x) \ <P C 
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= C I Det(x)). (44) etc. 

~ ~ A systematic way of obtaining separate contributions 

Here | Det(x)) = E^x) E A {x) \ Det) is a new Slater £ [n] , g) . g ted b the UCQA (gee alg0 m 

determinant composed ot the renormaiized orbitais, , n , , , , . „ 

where one can tind general analytic expressions lor the 



\cj> a {a;x))=E a (x)\(j> a {a)) [a=\,...,A), (45) 



corresponding correlated operators). In case of commuting 
operators f(a,/3) (e.g., for the central correlation factors 
/(a,/3) = l + h(\ r a ~rp |) depending only on the distance 
between particles) one can write (cf. Appendix A in |32|). 



J = cxp(^ ln[l + h{a, /3)]). (54) 

a</3 



\Deb{x)) = -j= e T V{\<f> pl (l;x))---\ct> PA (A ] x))} 

' V£S A — 

(46) 

In turn, such orbitais can be evaluated in a concise an- 
alytic form as initial ones are linear combinations of the Aft . ,i • i • ,i TTr , A1 ,, , , 
ttai r , ., , , a j • a \ After this, applying the UCOM procedure we get 
HOM orbitais (see Appendix A). 

Following ( 40 1 we arrive to /\ > \ n i £t/ q m r £.t / \p / \ 

B ' <9^(y) = [i + ^ T (a,/3)]{£ ; ( ! t (-y)-Ec e (y) 

F corr (x,y) = ($ corr (-x) | Eli-y^iy) \ $ corr (x)) = +^J(-»)^(y)}[l + fc(a,0)] 

= {Det(-x) | C ,t £ , J(-y)Sid/)C' | Det(x)}. (47) -El(-y)E a (y) - e\ {-y)\E p {y) (55) 



2 2 

FcorrK*) = exp(^) J P C (g,w) ) (56) 



Expressions ([26]) and (J29) with expectations F(v, s) and Along such a guideline we obtain putting in eq. |48| once 

F(w ', s '), which are determined by eq. ( 47 1 , are certain x — v and y = s by eq. ((28 I 
base for our calculations. 

2.2 Calculations with the Jastrow-type correlator 

I A 

We have seen how expectations ([22]) and ([25 1 with respect i*b(g, u) = -7 (Det(-v) | ^ e^^C | Det(v)) , (57) 

to the correlated WF (35 1 can be expressed through the a=i 

generating function 1 ,1 _ 11 1 u /ToTT 

6 6 and then x — v and y — s by eq. (31 1 



F corr {x,y) = \{Det{-x) \ Q corr {y) \ Det(x)), (48) 



2 2 



^ v , ™<^. v y F corr K, S ')=exp(^)iV c ^, V '), (58) 

Qcorriy) =&J2 El(~y)E a (y)C. (49) I ^ 

a— 1 

Since we are going to demonstrate (at least, qualitatively) (gg^ 
the CMM effects on the FFs and MDs against the SRCs when deriving these formulae, we have used the relation, 
inclusion ( 35 ) , let us employ, first of all, the Jastrow ansatz 

exp(-y*a Q f ) exp(-ya Q ) = e^ y y exp[-y*a Q 1 + ya a ], 
C = J = /(1j2)/(1,3).../(1, A) t j iig spec jg c realization of formula ([l8| for any c-vector y. 



x/(2, 3).../(2, A) Our consideration is simplified if Det(x) becomes in- 



xf(A-l,A). (50) 



dependent of the vector x, i.e. 



Then we have the decomposition \ Det(x)) =\ Det(0)) =\ SD) , (61) 

Qj(y) = J^Q^(y)J = Q^(y) + Q^(y) + where | SD) is an original Slater determinant (see below). 

_ Then 
+Q [A] (y), (51) 

1 A 

where (y) is an n-body operator so that ^c(<7, u) = ^c(<?, 0) = ^-(5D | ^ e^C | SD) 

a=l 

A (62) 

QW(y) = ^£;t(_ y ) jBa ( y ) ; (52) and 

a=1 A 

a N c (z,v') = N c (z,0) = ±(SD\&J2 e iz ^C | SD). 

Q [2] (y) = E^(y)' ( 53 ) 

a <( 3 ( 63 ) 
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In accordance with eqs. (22 1 and (24 1 the corresponding 

(64) 



and 



FF and MD can be written as 

F E ST(q) = F TB (q)F c (q 
with 



Fc(q) 



(SD | CV^C | SD) 



and 



(SD | CtC* | SD) 
</wd<> r^w j e-^ z N TB {z)N c (z) l lz 



(2tt) e 



with 



(5D | &e ix Pi6 | 5D) 
<SI> | CtC | 51?) 



The canonical TB factor 



Ftb(q) = exp( 



£^0> 
4A ' 



(65) 



(66) 



(67) 



(68) 



has appeared in formula (17) for the intrinsic operator 
F int {q), while 



N TB (z)=eM 2 -^) 



(69) 



is the own TB factor (see discussion in ref. [15]) for the 
intrinsic MD. Respectively, the function Fc(q) and the 
Fourier transform 



Vc(p) 



(2 



lyT / e-^N G (z)dz 



(70) 



determine the no CM corrected FF and MD with the cor- 



related g.s. (35 1 normalized to unity. 



To go on our exploration with Jastrow-type correla- 
tions, let us write down instead of eqs. (57 1 and (59 1 as in 
eq. ffi), 



Fj(q, v) = ( g , v ) + F ™ ( g , v ) + ... + F W (q, „) ( 7 1) 
and 

Nj(z, v') = (z, v 1 ) + [z t v ') + ... + nW ( Zj „') 

(72) 

to obtain with the help of the UCOM the following ex- 
pressions: 



a=l 



1 A 

F^{q,v) = -^Det{ r v)\Y t \P^P)-^ 



a</3 



:[e iq9 " +e iq ^]\Det(v)), (74) 



1 A 
N^(z,v') = j(Det(-v') | ^ I Det(v% (75) 

a=l 

N^{z,v')= l x (Det{-v')\ 



A 



J2{f(a,(3){e^+e^}f(a,(3) (76) 
-e !Z? « -e lz P?} | Det(v')), 



for central correlation factor f(a,(3) — f(\ r a — rg 
(a,p = l,...,A). 



2.3 Application to 4 Ue 



On the condition (61 1 the matrix elements ( 73 )— ( 76 1 are 



transformed into the corresponding expectations with re- 
spect to the | SD) . Such a situation is realized for the pure 
HOM (Is) 4 configuration occupied by the four nucleons in 
4 He. Indeed, it is the case, where the orbitals 

I 0a («)) =1 I X<rr(a)), 

is annuled with the operators a a (a = 1,...,4) so 
the renormalized orbitals (48) coincide with the initial 
| cj)a(ct)). Here x&t is the spin (isospin) part of the orbital 
(err = ++, H — , — h, )• In other words, the correspond- 
ing determinant (46 I does not depend on x, i.e., 



Det{x)) =| Det(0)) =| (Is) 4 ). 



(77) 



Taking into account the definitions ( 65 ) and ( 67 ) , the 
quantities in question can be represented as the ratios, 



and 



where 



F C (q) = Fj(q) 



N c {z) = Nj(z) 



AM 
Aj(0) 

Bjjz) 

Bj(oy 



i A 

F ll] (q,v) = - (Det(-v) | V e lq? ° | Det(v)), (73) and 



A/(q)H(ls) 4 |jV« ? iJ|(ls) 4 ) 

= A^(q)+A^(q) + ...+A^( q ) 



B, 7 (z) = ((Is) 4 | jV^J|(ls) 4 ) 

= B^(z) + BW(z) + ... + B^{z) 



(78) 



(79) 



(80) 



(81) 



so that Bj(0) = Aj(0). 

One should point out that we prefer to deal with finite 
decompositions (80 1 and (81 1 retaining for our approxi- 



mations only a few first terms of them. Effects of the ne- 
glected terms can be estimated (at least, for 4 He as in [43J) 
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by means of a direct computation without any decomposi- 
tion (see sec. 5). Of course, the numerator and denomina- 
tor in each ratio ( 78 1 and ( 79 1 should be equally truncated 



to meet the requirements FjJO) = 1 and Nj(0) — 1, which 
guarantee the correct normalization of DDs and MDs. In 
the context, we will recall many works based upon the 
so-called ?y-expansion (see paper [5] and refs. therein) of 
the inverse denominator Aj (0) in a series. In our opin- 
ion, such a procedure create some problem of convergence 
even for finite A. 
Thus we assume 



Aj( q ) = A^( q ) + A^( q ), 
Bj(z)=bW{z)+bW{z) 



with 



ip 2 s (r)e^dr, 



B^{z) = {e^) = J rt s ( P y zp d P , 



A-l 



(M(q)) 



(82) 
(83) 

(84) 
(85) 

(86) 



a</3 



and 



A-l 



(B 12 (z)). (87) 



a</3 



Here 



A a p{ q ) = exp[-iq(r a +rp)} 



x{[/>,/3)-l]e^«-^> +H.c.}, 
B afj (z) = exp[~iz(p a +p p )] 

x {f(a, p)e^ z{ P-'P^f(a, f3) - e¥*® a -v?) + H .c.}, 



(89) 



(a, = 1,..., A) 



since [f(a,/3),r a +r ] = [f(a,0),p a + Pp] = 0, and the 
symbol (...) is used to denote the expectation with respect 
to the determinant | (Is) 4 ) (generally a | SD)). In eqs. 
( 84 )— ( 85 1 fi s (r) (if ls (p)) is the Is orbital in coordinate 
(momentum) representation. For convenience, the general 
HOM orbitals are given in Appendix A. 



Further, calculations by formulae (B.3 1— ( B.8 ) with the 
HOM orbital (fi s and the correlation factor QB.12 ) are re- 
duced t o si mple quadratures. In particular, the approxi- 
mation (82 1 results in the FF, 



Fj(q) 



Ml) 

Moy 



(90) 



Aj{q) = aiexp(- 



462 



1 1 

"2exp(-^|) +a 3 exp(-^|, 



(91) 



with the coefficients 
ot\ = 1, a 2 = — 



6 



(l + 2y) 3 / 2: 
and the falloff parameters 



«3 



(1 + 40)3/2 



b i = r a = Po, h^hJ — 



l + 2y 



l + Ay 
l + 2y' 



bi < b 2 < b 3 



The DD associated with FF (90 1, i.e., its Fourier trans- 
form, can be represented as 



Pj{r) 



n-*/ 2 b\ 

Mo) 



[di exp(-6 2 r 2 ) + d 2 exp(-& 2 r 2 ) + d 3 cxp(-6 2 r 2 )] 



(92) 



di = 1, d 2 = - 



(1 + y) 3 / 2 ' 3 {l + 2y) 3 / 2 ' 



At the same tim e the approximation ( 83 ) gives rise to 
the MD (cf. eq. ((701), 



vj(p) 



1 



(2tt) £ 



e- ipz Nj(z)dz 



MO) 



1 p 2 1 p 2 1 p 2 

x ex p(- - ^) + fa exp(- — ^ ) + & exp(- - ^)] , 

(93) 



with 
and 



(l+3y)3/a 



[(l + 4y)(l + 2y)]3/ 2 ' 



7i = l, 72 = ———, 73 = 1 + 2?/. 
1 + 2y 

Henceforth we introduce the dimensionless parameter 

V = ( — ) 
r c 

The corresponding CM corrected quantities are deter- 
mined by 

Fj,EST(q) = F TB (q)Fj(q), (94) 
PJ,EST(r) = J e-^ r F JiEST ( q )dq, (95) 



so 



Vj ,est(p) = J e- ipz N TB (z)Nj(z)dz, 

^/ 2 b\ 

PJ,EST{r) = 

x [di exp(-6 2 r 2 ) + J 2 exp(-blr 2 ) + J 3 exp(-6 2 r 2 )] 



(96) 



(97) 
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with 



where 



h = 



VI- A- 1 ' 



h = 



and 



VJ,est(p) 



b* = 



^- 3 / 2 z>r 3 

4,(0) 



1 p 2 



1 p 2 - 1 p 2 - 

[/3iexp(- — + ^ 2 exp(- — -~) + /? 3 exp(- — —)] 
71 Ol 72 »i 73 of 



(98) 



with 



/ \ 3/2 / x 3/2 , v 3/2 

V7i/ V72/ V73/ 

where 

71 =1- A -1 , 72 = 72 -^ _1 , 73 = 73 -^~ 1 - 

For this consideration the mass number A = 4. But we 
preserve in these formulae the A~ l - dependence to indi- 
cate a distinct feature of the CMM correction. 

The approach developed here can be useful even if the 
ansatz (61 I does not work. In particular, the separate con- 



tributions to the matrix element ( 48 1 , that stem from de- 



composition (51 1 j can be expressed through certain s.p. 
overlap integrals with arbitrary orbitals <f> a . For example, 
we have for the one-body contribution (73 1, 



whence 
1 



F, 



tu) <(q,v) = <¥>!,(-«) I I ^»)K s (-u) I cp'MY 



(99) 

where | <p' ls (a;v)) = E a (v) \ <pi s {a)) (a = 1,2,3,4) the 
renormalized s.p. state and omitting the label a we denote 
I <P'is( v )) = E(v) I (pis) (cf. eq. (45 1). Analogously, one can 
get 



Flf s)i (q,v) = -( l p[ s (l;-v)\(^ s (2;-v)\ 
x If {1,2) - l][e lq91 +e lq92 } 
x|^i,(l;t;))|^(2;t>))^i,(-«)l^a(«)> 5 



(100) 



Let us stress once more that if the vector | tp ls ) is a lin- 
ear combination of the Cartesian states | n x n y n z ) the s.p. 
matrix elements involved are calculated using purely alge- 
braic means. 

In addition, we would like to show some results ob- 
tained with the Darmstadt (D) correlator, which is de- 
termined by eqs. ( 37 1— ( 39 ) . It is the case, where, e.g., 



instead of the operator Ai2(q) in expectation (86) one 
should write, 



«3(1>2). 



i lqf +H.c.} (101) 



For brevity, we introduce the CM coordinate R = i(ri + 
r 2 ) of particles 1 and 2 with their relative coordinate r = 
n — r 2 and momentum p = \{p\ — p 2 ). 

The hermitian generator used in |33| looks as 

3(1,2) = h-^lrp + pr^} (102) 
2 r r 

One expects the unitary operator c = exp[— ig(l, 2)] to 
shift the relative distance r between the particles via the 
position-dependent displacement s(r). A key point is to 
find an appropriate function s(r) such that c(l,2) could 
be tractable as a correlator in coordinate space. In the 
context, the authors of work have shown that 



FW(q,v) = I^ 0l (l;-«) I (cf> a2 (2;-v) \ ...{<j> aA {A--v) \ 

A 

xJ2e iq9a A A I ^ 0l (i;«)> I Kftv))- I <t>*AA;v)) 

a=l 

Remind that the antisymmetrizer 

A A = e vV 

ves A 

acts onto the subscripts of orbitals </> a (o:; x). It means that 
for (Is) 4 configuration 

F^ s)i (q,v) = ±(<t> ++ (l;-v) \ (<f>+-(2;-v) \ (</>- + (3; -v) \ 
4 

x (^__(4;-«) | $>^M 4 I cfi ++ (l;v)) | </>+-& v)) 

a=l 

x I ^_ + (3;v)} I ^ — (4;v)> 



f fl = 2(l,2)fc(l,2) = 



R + (r) 



r. 



(103) 



where the shift R+{r) — r characterizes some deviation of 
the transformed distance r g from the uncorrelated original 
r. 



The relationship (103) enables us to write 

AM?) 



ct(l,2)e^c(l,2) = ex Pl -i- 

2 r 



qr] 



(104) 



Substituting (1041 into eq. (101 1, we obtain with the 
(Is) 4 configuration, 

A [ $( q ) = \{A?M) = 3 eXp[ -j_ g2r ° ] C (g ), (105) 



C(q) = 



Stt 



.:>,.. i,„ir» r sin hl R +( r ) sin§qr 



r 2 dre 2P ° r { 



R+(r) 
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The property C(0) = provides the required value 
Fd(0) = 1 of the corresponding FF, 



Thus 

E(v) | 1pm) =| 1pm) + v m \ Is), (m = 1, 0,-1) (110) 



F D {q) 



A D {q) 



A^(q) +Ap(q) (106) Obviously, the second term in the r.h.s. of eq. (110 1 does 



A D (0) 

Furthermore, one can find the relation, 

l 

R + (r) = r + $(1,2; sd r )s(r) = r + J duexp(usd r )s(r) 

o 



As anticipated, f or a s mooth shift function s(r) small com- 
pared to r from (1071 it follows (cf. eq. (62) in |33]). 



R+(r) = r + s(r) + 



(108) 



One should note that the authors of [33] not indicating 
any model for s(r) have preferred to work with the cor- 
relation function R+(r) directly. Our calculation with a 
parameterized (sophisticated) form for i? + (r), taken from 
[33 , will be presented somewhere else. 



not contribute to the determinant | D(y )) that immedi- 
ately gives rise to (1091. 



As before, such an observation essentially simplifies our 
consideration since the matrix elements (73 1— (761 and so 



on are reduced to the expectations with respect to the 
customary shell determinant | ( ls) 4 ( l p) 12 ) . Owing to this, 
one can again employ formulae (B.3HB.8) to get the FFs, 



( 107 ) DDs, and MDs without any CMM correction 



Fj(q) 



AM 

Aj(oy 



(in) 



Aj{q) =ai(g)exp 



46 2 



+ a 2 (g)exp ( ) + a 3 (g)exp 



46| 



(112) 



2.4 Application to 16 

For another j-closed nucleus le O we will start with the 
fully occupied (ls) 4 (lp) 12 configuration which is built 
from the corresponding HOM orbitals in the Zs-coupling 
scheme (see Appendi x A ). Now, all we ne ed i s to show 
that the relevant SD (46 1 has the property (61 1. In other 



words, let us verify the relation 

| Det(v)) = Ei(v)...E 16 (v) | (ls) 4 (lp) 12 ) =| (ls) 4 (lp) 12 ) 

(109) 

for any vector v. 

Indeed, along with the evident equation 

E(v) | Is) = e" a | Is) =| Is) =| 000) =| 0) 
we find step by step, 

1 , , , 1 „+ i 



I Ipl) = 

I 1P-1> 
and 



100)- 



010) = (- 



I lpo) =| ooi) = 4 1 o), 



7i |100) -7i |010 ^ ( 7i 4 -7i^ 



e° a | lpl) =| lpl)+«+i | Is), 
e" a | lpQ) =| l P 0)+v | Is), 

e" a \lp-l) =| | Is), 

with the cyclic components 
1 



v± = =F—f=(v x ± iv y ), v 



Vz 



0), 



0), 



pj(r)= 7r -^{di(r)exp(-b 2 1 r 2 ) 



d 2 (r) exp(-6^ 2 ) - d 3 (r) exp(-b 2 2 r 2 )}, 

(113) 



7r- 3 / 2 6r 3 1 p 2 

vjfr) = a ,^ mm ex p( — v?)+ 

Aj{0) 7i bf 

1 p 2 1 p 2 

/3 2 (p)exp( to) +/3 3 (p) exp( 

72 of 73 bf 

vs. the CMM corrected ones, 

Fj.EST(q) = F TB (q)Fj(q), 

7r- 3 / 2 6i 3 - - 2 2 x 

PJ,EST{r) = A , Q . [d x (r)exp(-6i r ) 

+ d 2 (r) exp(-6~ 2 V) - d 3 ( r ) exp(-6~ 2 V)], 

7r- 3 /267 3 - p 2 
VJ,est(p) = . , n , [ffi(p)exp( 



4/(0) 



2 2 
P - p 

/3 2 (p) exp(- l2 - r ) +/3 3 (p) exp(- 7 ^-)] 



Of course, here we have the relevant TB factor, 



FrB(q) = exp( 



2 2 

g rg, 
64 ' 



(114) 

(115) 
(116) 

(117) 



(118) 



Analytic (in general, cumbersome) expressions for the 
polynomials a;(g), d 4 (r), d~i(r), ft(p) and /3j(p) (i = 1, 2, 3) 
can be obtained using formulae of Appendix B that results 
in (by taking, respectively, x — qjb\ and z — p/b\) 



a x (q) = 1 — - , a 2 (q) = 2 



iru(y) + K 12 {y)x 2 + ir 13 (y)x 4 



(l + 2y) 



3/2 
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03(g) 



7T2i(y) + n 2 2(y)x 2 + ^23{y)x 4 



where 



7rn(y) = -1 
112(2/) = 
^13(2/) 



(l + Ayf 2 

50 + 116y + 77?/ 2 
4(1 + 2y) 2 ' 

13 + 2hy + 22y 2 + 18y 3 
8(l + 2y)3 

_ 2 / 2 (-3-2 2/ + 5y 2 ) 



16(1 + 2y) 4 

7T2i(j/) = Tli(j/)| y ^ 2 » (* = X > 2 > 3 ) 



SO 



di(r) = 1 



/> 2 r 2 /3 
1 ( 13 *,2„2 



2 V 2 



da(r) 



ITT ^( 7 r 11 ( 2/ ) + 2. 12 ( 2/ )(3-2 & 2 r 2 )(^ 



+ 47r 13 ( y )(l5-206 2 r 2 +46^ 4 ) (|) ) 



(IT ^(. 21 (,) + 2 7 r 22( ,)(3-2^ 2 ) 
+ 47r 23 (2/) (15 ^206 2 r 2 +46^)gj ) . 
At the same time we find for the MD, 



1 z 2 
/MP) = 4 + Y , 



72 7 2 / 



&(p) 



»»i(i/) + 4(|-^)(»?22(y)/'y3) 



where 



^n(y) : 
7713(2/) 



(l + 4y) 3 / 2 (l + 2?/) 3 / 2 

3 3ly 2 + Uy + 18 
~2 (l + 2y) 2 ' 

13 + 69y + 92y 2 + 42y 3 
4(l + 2y)3 

3 y 2 (1 + 4y + 3y 2 ) 



mi{y) = 



(l + 2y)4 
3 9 + 44y + 62y 2 



2 (l + 4y) 2 
13 + 34y 
~8(l + 4y)' 



Th e cu toffs 61, 6 2 , b 3 , 71, 72 and 73 are determined as in 
eq.(pll and eq.d93l. 



The analytic expressions for the polynomials e?j (r) and 
(i=l,2,3) are obtained by following the recipes: 



d 2 (r) = 
d 3 (r) 



1 - 



1 - 



1 l + 2y 
A l + y 

1 l + 4y 



-3/2 



d 2 (r)| 6l: 



-3/2 



*bl,2' 



>i>1.3' 



A 1 + 2y 

h{p) = Pi 0)1 71 ->-7i> 

= 72~ 3/2 /3 2 (p)l7i,2^7i,2. = 73~ 3/2 ^3(p)| 7 i, E 



3— ^71,3 ' 



where 61, 6 2 , 63, 71, 7 2 and 73 are determined in the same 
way as in the case of He. 



3 Results and discussion 

The analytic expressions derived in sect. 4 for density and 
momentum distributions and their Fourier transforms are 
sufficiently general to be applied in different translation- 
ally invariant treatments with the S RCs i nclu ded. Our cal- 
culations carried out by formu l ae (|9 0|)-( 98 ) for the 4 He 



nucleus and by formulae ( 111 h(U18|for^°0 nucleus are 



displayed in figs. 1 — 5 together with available data. In 
these figures we distinguish two cases in which along with 
the model Jastrow correlations the CMM correction is ei- 
ther included or not. 

In order to calculate the charge FFs we have used the 
relation 

F CH (q) = F TB (q)F DF {q)F proton (q)F mt (q) 7 (119) 

where Fnp(q) = 1 — q 2 /2m 2 is the Darwin- Foldy correc- 
tion and Fp ro t n (q) is the finite proton size factor with the 
parametrization from |55| . 

The parameters r and r c (or, equivalently, y — 
( r o/ r c) 2 ) have been extracted from the data in fig. 1 for 
each nucleus via a least squares fit to the experimental 
Fch(q) '■ their best-fit values are ro = 1.163 fm and 
y = 3.120 (r c = 0.658 fm) for 4 He and r = 1.171 fm 
and y = 20.192 (r c = 0.261 fm) for le O. Being fixed in 
such a way, they remain unchanged for subsequent calcu- 
lations. Along with the best-fit solid curves we have drawn 
the corresponding dashed curves to demonstrate the CMM 
influence (sometimes considerable) on the distributions in 
question. As seen in fig. 1, the CMM-corrected calcula- 
tions reproduce the observed q-dependencies of the FFs, 
viz., the envelopes of diffraction maxima and the positions 
of diffraction minima. 

In order to evaluate validity of the approximation given 
by eqs. ( 82 )-( 83 1 we have calculated quantities Aj(q) and 
Bj(z) without any truncation of decompositions (80 1 and 
(|81|) . Comparison between the corresponding curves shows 
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Without CMM 



CMM 



r =1,163 




r =1,163 


y=3,120 




y=3,120 






x 2 =2,5 


dashed line 




solid line 




Without CMM 



r,=1,710 




r„=1,710 


y=20,192 




y=20,192 


z ! =4,212 




X ! =1,897 


dashed line 




solid line 




q (fm ) 



q (fm"; 



Figure 1. The charge form factor of the nuclei 4 He (on the left) and la O (on the right) : calculated with the Jastrow WF using 
the EST prescription (solid curves) and without the CMM correction (dashed curves); experimental points from |58| and [59] . 
respectively. Other clarifications are given in the text. 



Without CMM With CMM 



r„=1,163 




r=1,163 


y=3,120 




y=3,120 


dashed line 




solid line 




r(fm) 




Figure 2. The charge density of the nuclei 4 He and le O: calculated with the Jastrow WF using the EST prescription (solid 
curves) and without the CMM correction (dashed curves). In addition, the thick solid and dash-dotted curves show our exact 
(numerical) calculation for 4 He, respectively, with the EST prescription and without it; experimental points from |56j. Other 
clarifications are given in the text. 



that some qualitative changes of the r— and p— dependen- 
cies pj^srir) and 7?j,.est(p), which are determined, re- 
spectively, by ( 97 1 and ( 98 1 , can be by-products of the ap- 
proximation. In fact, considerable dips in the solid curves 
on the left panels of figs. 3 and 5 do not appear for ex- 
act calculations. At the point, one should note that the 
additional depression of pj,est(t) with respect to pj{r) 
at a moderate y— value (c/. the solid and dashed curves 
in fig. 3 for the alpha-particle in the range < r < 1) 
is obscured in the charge density. The latter, being de- 
fined as th e Fourier transform of the charge FF by for- 
mula (1191, is calculated via the convolution of pj^ESTif) 



with a smoothed charge distribution in the proton. More- 
over, it turns out that even with the lack (at smaller 
y— values) of the necessary property of pj t EST{r) to be 



positively definite the convolution results in a distribu- 
tion Pch{t) which has much in common with that shown 
by the solid curve in the left panel of fig. 2. Perhaps, in 
spite of similar observations many authors (see, e.g., Ta- 
ble I in [TSj with the parameters b = r$ = 1.1732/m 



and f3 



2.3127 fm 2 for 4 He that is equivalent to 



y=3.183) show only the charge densities of nuclei. Further, 
the exact distribution p < j X E < ST( r ) (^ ne thick solid curve in 
fig. 3) has a plateau in the vicinity of r = with a shallow 
dip. When increasing the y— values the pj(r) dependencies 
(both exact and approximate) become smoothly varying 
functions of the nucleon coordinate r. 

In addition, as seen from figs. 3 and 4, the CMM cor- 
rection diminishes the expected depression of the intrin- 
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Figure 3. The point-proton density of the nuclei 4 He and ie O : calculated by formulae (92 I (dashed curve) and (971 (solid 
curve) on the left and by formulae (1131 (dashed curve) and (1161 (solid curve) on the right. Distinctions between other curves 
are the same as in fig. 2. The normalization is J pj(r)dr = 1. 



0.16 




r(fm) p(fm" 1 ) 



Figure 4. The one-body density (on the left) and momentum distribution (on the right) of the alpha particle at different 
y-values and fixed ro = 1.163/m. As in fig. 3, curves on the left panel calculated by formulae (92 I and (971 whereas the right 
panel demonstrates the dependence 77 j,est(p) (eq. (98 1: solid curves) vs. rjj(p) (eq. (931: dashed curves). Distinctions between 
the thick solid and dash-dotted curves are the same as in fig. 3. The normalization is J rjj(p)dp = 1. 



sic DD pj(r) relative to PhomI 1 ") = lim^co pj (r) in its 
central region, i.e., increases the probability to find a nu- 
cleon in the 4 He interior. From the physical viewpoint 
such an extra increase is not something exclusive since the 
TI restoration means the introduction of nucleon-nucleon 
correlations ctS £l whole (including the short-range ones 
too). 

Going on our discussion of the interplay between the 
CM fixation and the phenomenological introduction of 
N — N repulsion in the nuclear wave function, we will 
note a simultaneous shrinking of the OBDD and OBMD 
(cf. the thick solid curves vs dash-dotted ones in figs. 3, 
4 and 5). Following [15 the term 'shrinking' implies that 



the EST prescription gives rise to increasing each of these 
densities in their central regions (respectively, < r < ro 
and < p < po — r^ ) compared to the nTI quantities. 
But unlike refs. [15] and |57|, where the effect has been 
confirmed within the HOM and its modification [S3] , the 
present observation is related to the exact numerical re- 
sults obtained beyond such simple models. In the context, 
note the relations under the strong inequality r c <C 
with 

pfwMir) = lim P j,EST(r) = f^ 3 exp(~r 2 /fg) 

y-s-oo 

VS. 

pHOM{r) - r^ 3 7T- 3 / 2 exp(-r 2 /r 2 ). 
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Figure 5. The momentum distributions of the nuclei 4 He and le O. Together with our calculations for the best DDs we have 
depicted the results from |60| (circles on the left) and |61j (the dash-double-dotted curve on the right). Difference between the 
thick solid and dash-dotted curves explained in fig. 3. 



and 



In its turn, the CMM correction contributes to 



Vmm(p) = lim v J,est (p) = r§7r 3/2 exp(V/Po) 



vs. 



lim r)j(p) = Vhom(p) 



-3/2 



exp(-p 2 /P 



Remind that here r = \/3/2r andpo = V^/2po ^ f^ 1 so 
we see one and the same renormalization of the parameters 
r and po m accordance with the conclusion from [15 that 
the so-called Tassie-Barker factors should be different for 
different distributions of particles in finite systems. 

Now, one can ask to what extent the mean square radii 
of these DDs are modified due to the CMM corrections and 
the SRCs effects. The analytical expressions of the FFs 
obtained here enable us to find an explicit dependence of 
the corresponding radius on parameters tq (po = Tq 1 ) and 
y. In this connection, let us recall that it can be found as 
coefficient of —q 2 /6 in the conventional expression F(q) = 
1 — lq 2 r 2 ms + • • • . In particular, we get 



Fj(q) 



i-iflV) J + . 



with 



r(0) 



= -6^ 

4/(0) 



(120) 



where Aj(q) is given by eq.(91l (eq.(112)) in case of 4 He 
( 16 0). Here A' 7 (0) = ^Aj(q) | g=0 . Doing so, one can 

(r 2 )j - (r 2 ) HO M, where 



evaluate the difference Aj 



IHOM 



3„2 (9 2- 
2 r \4 r 0, 



for He( O). For example, Aj = 
0.282 fm 2 at r = 1.163 fm and y = 3.120 in case of 4 He 
and Aj = 0.195 fm 2 at r = 1.710 fm and y = 20.192 in 
case of 16 O. It means that along with the aforementioned 
depression the SRCs inclusion results in broadening the 
OBDD. 



Fj,Esr{q) = 1 - ^q 2 (r 2 )j.EST 



with 

where ctb 
expression 



\r )j,est = ctb 

,.2 



(r 2 )j, 



{r 2 ) hom / A. These quantities enter the 



)CH 



cdw + (r ) P + (r )j. 



that determines the rms charge radii (/ 
tracted from 



EST 



2\l/2 
ICH 



to be ex- 



FcH(q) 



i - \q 2 {r 2 )cH + 

b 



Remind their experimental values: 1. 676(2. 730)/m for 
4 He ( 16 0), taken from [55]. One can verify that these val- 
ues are reproduced by our calculations with (r 2 )j t EST = 
2.000 (6.644)/m 2 for 4 iJe( 16 0)) ((r 2 )j is equal to 
2.31 fm 2 and 6.77 fm 2 , respectively). Note also that ac- 
cordingly the prescription [55] (r 2 ) p = 0.775 fm 2 . 

Finally, one should note that we do not attach great 
importance to a fair agreement of our calculations with 
the data in figs. 1-2 and not too good one in fig. 5. In fact, 
as mentioned in sect.l, the I A, in which the charge FF is 
determined by formula (119 1, is insufficient (see, e.g., [11J) 



to give an adequate treatment of the elastic electron scat- 
tering off nuclei with the q-increasing when MEC effects 
become more and more important. In addition, one has to 
account for the hi ghe r -ord er contributions to the decom- 
positions by eqs. (80l-(81l. Once more it illustrates fig. 6, 



where we can see a considerable shift of the first diffrac- 
tion minimum towards the larger q- values. Of course, the 
shift may be compensated by modifying the values of the 
parameters involved. 
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q (fm" ; 



Figure 6. The exact point-like FF (thick solid line) for 4 He 
vs its approximation (solid line) by eq.(94l 



Along with the pronounced flattening of the thick solid 
curve in fig. 5 in the vicinity of the p — 2fm~ 1 it means 
that every time higher-order correlations effects should be 
investigated separately (cf. similar results obtained in [43, 
for A He without any CM corrections). In the context, the 
large difference between the Argonne [61J and our calcula- 
tions in fig. 5 at the p - values > 2 fm~ l can be explained 
to great extent by the inclusion of many-nucleon corre- 
lations in the former. Their role becomes stronger with 
increasing the mass number. In addition, being aware of 
the necessity [61 J of introducing noncentral correlations 
(see also [2]), we note that our method of restoring of the 
TI may be helpful for such complex numerical calculations 
as well. 



4 Summary 

We have shown how the approach developed in |15| when 
studying the one-body and two-body density matrices of 
finite nuclei can be realized beyond the independent par- 
ticle shell model. The appropriate treatment of the CMM 
is combined with the inclusion of the SRCs in the nu- 
clear WF, e.g., regarding either the Jastrow ansatz or 
the UCOA. In our translationally invariant calculations 
the OBDD and OBMD are expectation values of the A- 
particle multiplicative operators which are dependent on 
the relative coordinates and momenta (the Jacobi vari- 
ables) and sandwiched between intrinsic nuclear ground 
states. 

An algebraic procedure proposed earlier helps us to 
avoid a cumbersome integration and see certain links 
between the distributions in question being expressed 
through one and the same generating function. In the 
course of the procedure the so-called Tassie-Barker factors 
stem directly from the intrinsic operators (not the WFs). 
One can stress that these factors being different, unlike 
other works (see, e.g., [TS] and [5]), for the DD and MD 
occur by reflecting the translationally invariant structure 



of the corresponding intrinsic operators. Each of them is a 
Gaussian whose behavior in the space of variables is gov- 
erned by the size parameter ro(or its reciprocal po) an d the 
particle number A for a given finite system (nucleus) , but 
it does not depend upon the choice of the g.s. WF. The 
latter can be a simple Slater determinant, include SRCs 
or not, be CMM-corrected or not, etc. 

The use of the Cartesian or boson representation, in 
which the Jacobi variables are linear combinations of the 
creation d and destruction a operators for oscillator 
quanta, has allowed us to simplify the calculations for the 
closed shell nuclei A He and le O. Certainly, the underly- 
ing idea based upon the normal ordering of the operators 
that meet the Bosc commutation rules may be helpful in 
case of other closed and open shell nuclei. The analytic 
expressions for the intrinsic densities, form factors and 
momentum distributions derived in sect. 2 with the Jas- 
trow correlators are convenient in getting a deeper under- 
standing of some nuclear properties. In particular, after 
restoring the TI on the SRCs background we have both in 
pj{r) and rjj(j>) their shrinking at enough large values of 

the ratio y — (jr^J ■ 

Finally, regarding prospects of our approach in describ- 
ing the interplay between the CMM and the SRC effects 
we mean, first of all, its application for calculations of the 
two-body momentum distributions in such reactions as 
4 He(e,e'NN)X and 16 0(e, e'NN)X (cf.the correspond- 
ing qualitative findings in [15). Our work in the subfield 
is in progress. 



A A key point of calculations beyond HOM 



The algebraic technique, shown in sects. 2 and 3, can 
be also helpful in calculating the expectations by eqs. (22 1 



and (25) (or something like this) with WF <P that is either 
a linear superposition of SDs or a SD which is composed 
of (HF) or other model orbitals expanded in the HOM s.p. 
states. We find such expansions, e.g., for HF solutions [53_ 
and an effective inclusion [S3] of short-range repulsion be- 
tween nucleons (in both cases in spherical representation) . 

By definition, the normalized RKB-orbital (for a Is 4 
configuration in 4 £fe nucleus) is 



\<t>s) = 



1 



Wis) + P\<f>2s)) 



(A.l) 



with an adjustable parameter (3 . In this connection, let 
us recall the well-known expressions for the HO orbitals 
| nlm) that are specified by the principal (spectroscopic), 
orbital angular momentum and its projection quantum 
numbers n, I and m. One has in coordinate space 



Pnlm(r) = (r \ nlm) = R n i(r)Yi T 



(A.2) 



Rni(r) = C nt r 



3/2 ( r_ 
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<T> ( I //./ +- l;^) exp (-~s i 



2'ri 



C, 



V2 



nl 



r(i + l) 

while in momentum space, 



r(i- 



r(n) 



1/2 



<Pnlmip) = {P\ nlm) = Rnl(p)Yi m 



(A.3) 



3/2 ( P_ 

Po 

2 



where following [62] <P(a, c; x) is the confluent function. By 
passing, remind also the link with the associated Laguerre 
polynomials, 



B Relevant calculations 

The expectations of interest can be expressed in terms 
of these orbitals (in general, the s.p. orbitals |A) occupied 
in the g.s.) in different ways. For example, using the for- 
malism of secondary quantization, one has 

A®{q) = 23^ CT r£ Al , A2 eF<AiA 2 | A 12 (q) 

x | AiA 2 -A 2 Ai), (B.l) 

and 

x I AiA 2 - A 2 Ai), (B.2) 
where F means the Fermi sea, so 

AM(q) = A [ l( q )-AM c ( q ), (B.3) 

A 4iM = 4 E {<Px 1 <PX,\A 12 {q)\ V>XlV > Xa ),(BA) 



eY 2 (*) 



, <P ( l-n,l+ -;x ) , n= 1,2,. 



r(z , 2 

In turn, we find in the Cartesian representation 

\<P2s) = ^ \ n xn y n z )(n x n y n z \ip 2 s) (A.4) 

n x -\-ny-\-n s — 2 

one can show (cf.,[49 ) 

1^2.) = -i (1200) + |020) + |002)) (A.5) 



2 

41(9) = J E (VA^Aa I iiiafa) I ^^.(B.S) 



Ai,A 2 SF 



and analogously 



(B.6) 



B M = ^ E (<PXiVX 2 \B 12 (z)\ VXlV > X2 ),(B.7) 



Ai,A 2 £F 



or taking in account eq.( 16 1 

1 

\f2s) : 



— a) ■ d f |000), 
V6 



(A.6) 



i.e., for the RKB-orbital, 

\4> s ) = [1 + /^-^[l - G8/V6) a f a f ] | 0). (A7) 

Substituting (A. 7) into (A.6) (when calculating the ratio 
A IPM (q)/A IPM (0), the normalization factor [1 + /3 2 ]- 1 / 2 
can be omitted) we find 

exp ( X ■ a) | 0.) = [1 - 08/V6)(a f + *)(& f + x)] I 0) 

(A8) 

for any complex vector \. 

Now, after modest effort we obtain 

{<j> 8 | exp (-x * ■ & f ) exp (% • a) | S ) = 



= l + /3 2 - 3 /3 2 X*X- 
/3 2 

Ix'x' + nR y (x'x'Kn) (A9) 



B™ c (z) = j E (^a^a 2 |S 12 (z)|^^ Ai ).(B.8) 
Ai,A 2 eF 

We take the Zs-coupling scheme with the orbitals 

|A) = \<p X ) War) ■ (B.9) 

Accordingly eqs. ( [88p9| 

= ft f (l, 2) [e"? ?1 + e*^>] fc(l, 2) 
+#(1,2) [e l9?1 + e 4q?2 ] 
+ [e^ 1 +e"^ 2 ] ft(l,2), 



(B.10) 



B 12 («) = #(1,2) [e lz Pi + e**^] fc(l,2) 
+#(1,2) [e lz Pi +e 4Z ^] 
+ [e izg i +e lz Pi] h(l,2), 



(B.ll) 



once f(a,f3) = 1 + h(a,/3) (a,0 = I,..., A). In this 
work calculations have been carried out with the state- 
independent correlator 



h{a,/3) = h(\r a -rp\) = - exp 



(r a - rpY 



(B.12) 
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where r r is a correlation radius. 



Further, putting in the relation (cf. eq. 60 I 



exp 



— it a + ua 



= e-i u "exp 



— it a' 



the vector it equal first to 



and second to 



. r 
u = i—^q 

V2 



Po 

it = — =z 
V2 



exp (ua) 
(B.13) 
(B.14) 



we split expo nents exp (iq r) and exp (izp) , respectively, 
in eqs. (B.10) and (B.ll) into such a norma lly o r dered 
form. Then, when evaluating the sums in eqs. (B.4 1— ( B.5 1 



((B.7HB.8)), it suffices to consider the matrix elements: 
xe™ 1 | tpx^), (B.15) 

xe"" 1 | ip\ 2 <fx 1 )(k = 1,2), (B.16) 



we have employed the property ( 41 1 and introduced the 
operators 



H {1 \r;u) = h 



r H —it 

V2 



(B.17) 



and 



H&(r;u) = H ( - 1 \r;-u*)H^(r;u) (B.18) 



dependent on the distance r = n — r 2 between the nucle- 
ons. Obviously, the superscript k in (r; it) labels the 
order in the correlations involved. 



Using the definition (A. 2 1 and the transformation 



(110 1 the contributions of interest can be represented as 



„2 I -2 



xe ~%~H w (r;u) = 7 (fc) (V), (B.19) 

r 2 +r 2 

^ E m m[%. u (u) = r " 6 J dn / dr 2 e~ "VifW(r; «) 

x - «*n + ^ ri r 7 2 ) = J(fc) («) 

(B.20) 



tt 3 M«(m) = r - 6 / dn / dr 2 e 



2 , 2 
r l ~ ' 2 



2„,2 ^,2 i Jl 



H^(r;u) x 4{^ _ [i^^ 
r o 4r o 



2-v/2r§ 
(fe = l,2) 



[(rl + rl) (n - r 2 ) - r?n + r 2 2 r 2 ] } = («), 

(B.21) 



while 



K 3 Y,M { i k sa P m(u) = ^ 6 J dn J dr 2 



ri it 

— + -p 
r V2 



(B.22) 



% " 2 o H^(r;u)=I^(u), 



^E^W>) = 2r o~ 6 / dn J dr 2 



n it 

r \/2. 



x — e "o H {k \r-u) = J {k) (u), 
ro 

(B.23) 



^ 3 E M[% 1 . lpm2 {u)=Ar^ J dn J dr 2 e~ 

77117712 

xg( fc )(r; M ){ (T ^- (rir2 ^ V2) 

(rir 2 )(MT- 2 ) _ (ur 2 )(it*r 2 ) 1 = - (fe) 

3 _ 2r, 2 1 W ' 



\/2r3 



o 



k = 1,2. 



(B.24) 

Substituting expressions (B.17 1— ( B.18 1 into these equa- 
tions, we find with the correlator ( |B .12 ), 



(:2 



1 »« a /W(tt) = -/dn/dr 2 e-'"?- r ?e-2"" 2 



(rl + r|) , 



(B.25) 



2 J (1) (?i) = - J dn J dr 2 e- r *- r *e- vr \-^ vur 



x (r 2 + r 2 - it*it + 



it — it 



r), 



e 2! 



(B.26) 

! PW(«) =- /" dn J dr 2 e- r i- r h- yr2 e-^ yur 



if 2 2 r? -{- Vn t 

x 4{rfr^ - - " it*it 



it — it 



2^ 



[(»•? +rl)r-rln+rtn]}, 



(B.27) 



and 



C2 



-2 / dn/ dne^-^V^V^ 



v-2 



n, 



(B.28) 



e s»« 2 jW(u) = -2 / dn/ dr 2 e- r i- r 'e- wr2 e- 



n 



r 2 , 



(B.29) 



•>" 2 P«(V) = -4/ dn / dne-^-^e-^e-A^ 

[ (rir 2 ) 2 _ (rir 2 )(M'*7-2) 



+ 



(rir2)(^2) (ur2)(w*^2) 



2r 2 



(B.30) 
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It is readily seen that the corresponding counter- 
parts of the second order, multiplied by the same factor 
— exp \—\y ( it* 2 + u 2 )] , can b e obtai ned from the inte- 
grals QB.25| 1-( |R27) and ( |B.28| -( pHo} by doing in their 
integrands the two independent changes: y —¥ 2y and 
yu — > y (u — u*). In turn, these integrals may be calcu- 
lated by addressing an auxiliary integral 

I (u; a, bx,b 2 ) =fdrij olr 2 e~< r ' +r -) 

e -yr 2 e -V2yur exp + 

in the vicinity of the parameter values: a = 1, bi = b 2 = 0. 
Indeed, we have 



/ (u; a, bi, b 



x exp 



[a(a + 2y)f 2 

B 2 (b - V2yu) 
8a 2a + Ay 



(B.31) 



where B = b\ + b 2 and b = i (£>i — b 2 ) , and after evi- 
dent differentiating (for instance, using ana lytic means of 
Mathematica) we get formulae (eq lll)-(eq 1181. 
It yields 



and 



41(9) 



M2(u)+2ReM£l(u) 



(B.32) 



M^ c (u) + 2ReMgl(u) 



(B.33) 



and 



M dir( U ) = E 



(B.34) 



(B.35) 



Now we will separate out the purely Is subshell, mixed 
Is — lp and purely lp subshell contributions assuming 

<'(«) = M^{u) + M%l(u) + AfW(«) (B.36) 

and 

€(«) = ^f{u) + + MW(«), [h = 1, 2) 

(B.37) 

with 



MW(u) = Mi s fc )( M ) = (lsls | ffW(f;«) | lsls), (B.38) 
^itW = E [ M Spm(«) + M Sl;l S («)] . (B-39) 



mim2 

and analogously for the bar quantities. 



(B.40) 
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